Condensates of p-wave pairs are exact solutions for rotating two-component Bose gases 
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We derive exact analytical results for the wave functions and energies of harmonically trapped 
two-component Bose-Einstein condensates with weakly repulsive interactions under rotation. The 
isospin symmetric wave functions are universal and do not depend on the matrix elements of the 
two-body interaction. The comparison with the results from numerical diagonalization shows that 
the ground state and low-lying excitations consists of condensates of p-wave pairs for repulsive 
contact interactions, Coulomb interactions, and the repulsive interactions between aligned dipoles. 
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Exact analytical solutions for interacting quantum 
many-body systems are very rare [l|, |2|. However, they 
are of tremendous interest since they may provide us 
with further insight into the correlations of quantum- 
mechanical many-body systems. Even rarer are cases 
where the exactly solvable quantum many-body system 
can also potentially be realized experimentally. A famous 
example is the celebrated Laughlin state (Sj and its gener- 
alizations [J] that describe the two-dimensional electron 
gas in the quantum Hall regime. In recent years, the 
advances with ultra-cold atomic quantum gases have sig- 
nificantly broadened the range of experimentally acces- 
sible many-body systems that are also solvable exactly. 
An example is the interacting Bose gas in one dimension. 
Here, the Bethe ansatz offers an analytic solution to the 
Lieb-Liniger model [5| , and the fermionization of bosons 
in the Tonks-Girardeau [6| limit has been observed [7|,|8J. 

Exact analytical solutions exist also for single- 
component dilute and weakly interacting Bose-Einstein 
condensates under rotation J9l4l7| . For a large class of re- 



pulsive two-body interactions, the exact ground state of 
iV bosons at angular momentum L results from project- 
ing the unique state where L particles carry one unit of 
angular momentum onto the subs pace w ith zero angular 
momentum for the center of mass 



IJ-lia. In this Letter, 



we generalize this exact solution to the interesting case of 
a two-component Bose gas [18, 19,|, and arrive at a very 
appealing result: As in the single-component case, the 
exact solutions at angular momentum L are states where 
L bosons carry one unit of angular momentum each, and 
as before one has to project out excitations of the center 
of mass. In the two-component case, however, isospin en- 
ters as a good quantum number, and eigenstates can be 
labeled by the number of isospin-singlet p-wave pairs that 
enter the wave function. Isospin was introduced in nu- 
clear physics by Heisenberg 20] and - like ordinary spin 
- is based on SU(2) symmetry of two-component sys- 
tems. The comparison with numerical results shows that 



for repulsive contact interactions. Coulomb forces, and 
repulsive forces between aligned dipoles [2l|, the ground 
state contains a maximum number of isospin-singlet p- 
wave pairs. This Letter also explains the recent study by 
Bargi et al. [19j, who found by numerical diagonalization 
that the yrast energy of two-component Bose gases is a 
simple function of angular momentum. 

We consider a harmonically trapped two-component 
dilute gas of N bosons of a first species I, and M bosons 
of a second species f- We assume that the interactions 
are perturbatively weak and of equal strength between 
intra-species and inter-species. This special case of equal 
scattering lengths is approximately realized in gases of 
^^Rb and gases of ^-^Na [22|- Let us list the conserved 
quantities involved in this problem. Besides the isospin 
T, the total angular momentum L, the angular momen- 
tum of the center of mass Lc, the total number of bosons 
A = M + N and the isospin projection T^ = {M - N)/2 
are conserved. There is no simple basis that reflects these 
symmetries simultaneously. In second quantization the 
conservation of L, A, T^, and T is straight forward, while 
the conservation of L, M, N, and Lc is most easily ex- 
pressed in the configuration space within first quantiza- 
tion. Wc will employ both pictures in what follows. 

Single-particle states (t>nim of the spherical harmonic 
oscillator have energies Eni — haj{2n + I + 3/2). We 
are interested in low-energetic states at high angular 
momentum. For weakly perturbative interactions, only 
single-particle states with no radial excitation (n = 0) 
can contribute, and we can limit ourselves to maximally 
aligned states with m = I. Thus, single-particle states 
with single-particle angular momentum I are 4>oiiiz) = 
z'exp(-|zp/2) and (j>ou{w) = (w' exp (-|w;p/2) for the 
bosons of the species | and f, respectively. We denote 
the coordinates of bosons belonging to the ], species by 
Zj = Xj + iyj {j — 1 , . . . , N) and and by Wj = Uj + ivj , 
(j = 1, . . . , M) for the f species. Here, x and y {u and v) 
are the two Cartesian coordinates perpendicular to the 



axis of rotation for the J, (f) species. In what foUows, we 
omit the ubiquitous Gaussians from the single-particle 
wave functions. Let 6], and SL create a boson in the 
state corresponding to z' and w\ respectively. The cor- 
responding annihilation operators are 64 and 6;f, and the 
creation and annihilation operators fulfill the canonical 
commutation relations for bosons. 

In the case of a single species of bosons under rotation, 
the ground state at angular momentum L consists essen- 
tially of L bosons carrying one unit of angular momentum 
each while the remaining N — L bosons carry no angu- 
lar momentum 12] ■ Small modifications of this picture 
are due to the preservation of the angular momentum 
of the center of mass. Remarkably, the two-component 
case is somewhat similar, and Bargi et al. 19] numer- 
ically found that the ground state consists entirely of 
single-particle states with angular momenta I = and 
1 = 1. For L < N < M, there are L + 1 such states 
(labeled by the number of bosons of the species I that 
are in the single-particle state with I = 1). What is 
the ground state for repulsive interactions in the space 
spanned by these states? To address this question for a 
two-component system, we recall the essence of Hund's 
rule: For repulsive interactions, the interaction energy 
is minimized for wave functions that are antisymmetric 
in position space. In fermionic electron systems such as 
atoms and quantum dots, this leads to a symmetric spin 
wave function. The same reasoning applied to the present 
case of a two-species Bose gas would require the isospin 
wave function to be antisymmetric, too, thus making the 
total wave function symmetric under particle exchange. 
The operator 
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(1) 



creates an isospin-singlet p-wave pair (i.e. L = 1 and 
T = Tz =0). This pair is antisymmetric in position 
space and antisymmetric in isospin space and thus totally 
symmetric under exchange. Thus, the interaction energy 
is minimized for condensates of isospin-singlet pairs. In 
general, we have N ^ M, and the ground-state wave 
function will also consist of unpaired bosons. The states 



N- 



A-L 



T-) [bU) (&lt) [B^) |0) (2) 

have angular momentum L, isospin T = A/2 — r, isospin 
projection T^ = {M—N)/2, and consist entirely of single- 
particle states with angular momenta I = and / = 1. 
Here |0) denotes the vacuum, and r = 0, 1, . . . min {L, N) 
is the number of isospin-singlet pairs. The isospin oper- 
ators are f^ = Ei=o(^iV't ^ ^u^ii)/2, T^ = T,Zo b^ht, 
and f2 ^T^fL+TziTz + l). The eigenvalues of f^ and 
T^ are denotes as Tz and T{T + 1), respectively. Let 
us understand the state ([2]) in detail starting from the 
right. The application of the pair operators i?^ to the vac- 
uum yields a state of 2r bosons with quantum numbers 



L = T,T = 0,Tz = 0. The operators &L yields the angu- 
lar momentum L we seek, increase the number of bosons 
to i + T, while keeping isospin T = Tz = {L — t)/2 a 
good quantum number. The application of the operators 
&L increase the number of bosons to A and keeps isospin 
a good quantum number. Finally, the desired particle 
numbers M and N (i.e. Tz = (Af — N)/2) results from 
the T- operators. For repulsive interactions the ground 
state consists of the maximum number of isospin-singlet 
pairs (i.e. r = min(L,iV)). Thus, the observation by 
Bargi et ai, together with an adaptation of Hund's rule 
for bosons and isospin symmetry leads to eigenstates ([2]) 
consisting of condensates of isospin-singlet p-wave pairs. 
Note that these arguments are independent of the details 
of the repulsive interaction. These are the main result 
of the present Letter. As in the single-component case, 
minor modifications of this picture are due to the con- 
served angular momentum of the center of mass. Let us 
contrast our results to BCS pairing in Fermi systems. In 
BCS theory, a weakly attractive interaction leads to the 
formation of Cooper pairs (i.e. pairs of fermions in time- 
reversed orbits). The resulting BCS ground state is a 
condensate of spin-singlet s-wave pairs, where the sym- 
metric configuration-space wave function optimizes the 
interaction energy. In our case, we deal with bosons and 
with repulsive interactions, and this modifies the picture 
accordingly. 

Let us compute the energies of the stat es pi ). We 
generalize and extend the results of Ref. [1J| to the 
case of two-component Bose gases and sketch the main 
steps. The Hilbert space T-L)^ of N identical bosons l 
at total angular momentum L is spanned by products 
GAi iz)ex^ (z)- ■■■■ eAfc (z) (with Ai + A2 + . . . + A^ == L) of 
elementary symmetric polynomials 



e\{z) = e\{zi,...,ZN) = 



Zj, Zj, 



l<ii<...<ix<N 



(3) 

Note that e\{z) carries A units of angular momentum. 
For two-component mixtures of N and M identical 
bosons with N < M, the Hilbert space at total angu- 
lar momentum L (with L < M) is the sum 



min(L,Af) 
A=0 



(N) 
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L-A 



(4) 



of direct products of the Hilbert spaces of each compo- 
nent. Products of elementary symmetric polynomials are 
linearly independent and form a basis. In the absence of 
interactions, all states in the Hilbert space are degener- 
ate. Perturbatively weak interactions will lift this degen- 
eracy. 

The two-body interaction for a two-component mix- 
ture of Bose gases can be written as 



^ = y^ VmV,r, = y^ Vm{Ar. 



/ ^ ^ra ^ in 



m>0 
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c„ 
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Here, Vm is a matrix element and 
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>, (z,-z,r(5..-5.^r, 


(6) 


l<i<j<N 




>; {w,^w,r\d^^-d^j^ , 


(7) 


l<i<j<M 




>; {z.-w,r{d.^^d^^r 


(8) 



B„i — 



l<i<N 
l<j<M 

arc the interactions between bosons of species ^, between 
bosons of species fj and the inter-species interaction, re- 
spectively. By construction, the interaction preserves an- 
gular momentum L (i.e. the degree of the monomial wave 
function it is acting on), is of two-body nature, and - due 
to its translationally invariant form (only differences of 
coordinates and derivatives appear) - preserves the angu- 
lar momentum Lc of the center of mass. Furthermore, the 
interaction ^ is invariant under the exchange of zi -o- Wk 
and therefore preserves isospin. For the zero-ranged con- 



tact interaction, we have Vm = {—l/2)"^/m\ [14 1 



The action of the operators ([6]) on elementary symmet- 
ric polynomials is particularly simple [IJ] 

Amex{z) = Bme^{w) ■ 
Crne\iz)e^{w) ■ 

AmBiiz) = B„iei{w) = C'raR ■ 

Here, R= ^ (ei(-z) + ei(w)) denotes the center of mass. 
Equations ^ show that only the terms < m < 2 of the 
Hamiltonian ^ are of interest when acting on products 
ex{z)ei_i,{w). For the operators Vq and Vi we find 

% = A{A -l)/2, V,=A{L-L,) , 






for m > 3 







for m > 3 







for m > 1 


(9) 
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L = ^Zidz 

i=l 
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Wjdnij , 
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N M 

^ ^ ZiWkd^-du 

.i,j = l k,l = l 



(10) 



Only V2 is truly non-trivial, and 

V2ex{z)e^,{w) = 2{XN + /lAf + 2Xfi)ex{z)ef,{w) 
+ 2{N-X + l){^i + l)ex^l{z)e^+l{w) 
+ 2{M-fi + l){X+l)ex+i{z)e^,-i{w) , 

- 2A{N - X + l)ex-i{z)e,,{w)R 

- 2A{M - fi + l)ex{z)ef,^i{w)R . (11) 

Equation pT|) shows that the set 

M = {R''ex{z)eL-x-n{w)} with (12) 

< A < min(i,iV) , <n< L- X 

spans a subspace in Hilbert space at angular momen- 
tum L that is left invariant by ¥2- This subspace con- 
tains the states ex{z)eL-x{z) which are linear combina- 
tions of the eigenstates ([2]). The states ex{z)eL-x{z), 



are, however, not eigenstates of the center-of-mass mo- 
mentum. Let Pq be the projector onto wave functions 
with zero angular momentum of the center of mass, i.e. 
PqiP{z, w) = -ipi^ — RjW — R) for any wave function 
ip{z,w). The wave functions Poex{z)eL-xiz) are in the 
subspace spanned by M. [1J|. Thus, the states PqIXt) 
with \xt) with from Eq. ([2]) are eigenstates of the Hamil- 
tonian ([5]). Note that these states do not depend on the 
matrix elements Vm of the interaction. 

To compute the corresponding eigenvalues we make the 
ansatz PoipL,n for the eigenfunction with 



min(L,Af) 

•ipL,n = ^ c^^'ex{z)eL-x{w) 

A=0 



(13) 



Here, n is an additional label that distinguishes between 
min(i, iV) -|- 1 different wave functions. Our results below 
suggest that n is the number of isospin-singlet pairs. The 
eigenvalue equation VPqiPl.ii = E„PoipL,n requires the 
coefficients c^' to fulfill 

= {LM + X{2L + N - M - 2X) - En) c^x^ (14) 
+ {N-X){L- X)c^;;_l^ + X{M -L + A)4"Ji • 

Here, e enters the energy eigenvalue 



En = A{A - 1) Y + ALvi + 2v2Sr, 



(15) 



Note that the eigenvalue problem (|T4l) does not depend 
on the matrix elements of the two-body interaction. For 
the solution of the eigenvalue problem, we make the 
ansatz 



.(") 



k=0 



(16) 



Here, we concealed the fact that the coefficients (3^ also 
depend on n. We insert the ansatz ([T5| into Eq. p^ 
and compare the coefficients of A™, m = 0, ...n -|- 2. 
This yields 



AL ~ n{A + 1 - n) , 



(17) 



which enters the energy P^ . The coefficients /3k, k < n 
are recursively defined in terms of /3„ (which sets the nor- 
malization). The quantum number n acquires the values 
n = 0,1,2, .. .inm{L,N), and the lowest energy is ob- 
tained for n = min(_L, N). 

Figure [1] shows the energy spectrum of a two- 
component Bose-Einstein condensate with iV = 4 and 
M — 8 particles per species, respectively, as a func- 
tion of the angular momentum L. The broken lines 
connect states with energies i?„ from Eq. P^ for fixed 
n = 0, 1, 2, . . . (from top to bottom). The ground state 



has n — mm{L,N). Note that c 



(0) 



1 solves the 



eigenvalue problem (fTil) and yields the state ipL.o 



cl^zi, . . . , zn,wi, . . . , wm)- This state is totally symmet- 
ric under the exchange of any particles, has maximum 
isospin T = A/2 and contains no isospin-singlet p-wave 
pairs. For L < N, the ground state p^ with n = L has 
coefficients c^Vcq^^ = (-l)^(M-L + A)!(iV- A)!/[(Af- 
L)!iV!] and can be rewritten as S']^^^-^(zfe—u)fc). Here the 
symmetrization operator S ensures the symmetry under 
exchange of bosons of each species, and the antisymmetry 
between the two species in position space is evident. We 
thus believe that the label n of the energies P^ has to 
be identified with the number r of isospin-singlet p-wave 
pairs of the eigenstates -PqIxt) with \xt) from Eq. ([2]). 



ing two-component Bose gases with weak repulsive inter- 
actions are condensates of isospin-singlet p-wave pairs, 
and we derived analytical expressions for the energies and 
the corresponding wave functions. The wave functions 
are universal as they do not depend on the details of the 
two-body interaction. Numerical computations demon- 
strate that these eigenstates are the ground states and 
some of the low-lying excitations for the contact interac- 
tion, the Coulomb interaction, and repulsive interactions 
between aligned dipoles. 




FIG. 1: (Color online) Spectrum of a two-component Bose 
gas with N — 4 and M = 8 bosons per species, respectively, 
as a function of the angular momentum L for the contact 
interaction. The broken lines connect states with energies 
En from Eq. (fT5)l for fixed n = 0,1,2, ... N (from top to 
bottom). The solid (red) line connects states with energies 
El for L<N. 

The reasoning that led to the isospin-singlet p-wave 
condensates ^ was based on general arguments re- 
garding repulsive interactions in SU(2)-symmetric two- 
component systems. To check our arguments, we per- 
formed numerical computations for Coulomb interactions 
and repulsive interactions between aligned dipoles. For 
the Coulomb interaction, the relevant matrix elements 
are vq = \/7r/2, vi = — wo/4, and V2 = Svn/di, 



re- 
spectively, and we refer the reader to Ref. [16| for the 
analytical derivation. The comparison with numerical 
results shows that the exact results (fTS]) again describe 
the ground states and low- lying excitations. Finally, we 
consider repulsive interactions between dipoles aligned 
perpendicular to the trap plane. The comparison of the 
numerical spectra and the analytical results (jl5|) yields 
vo « 6.868, vi sa —3.188, and V2 = 0.755, respectively. 
Again, the ground state and low-lying excitations are de- 
scribed by the analytical results. 

In summary, we showed that low- lying states of rotat- 
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